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Galaxy clusters furnish extremely rich information on the contents and structure of our uni- 
verse. The potential of galaxy cluster studies to constrain dark energy, for example, motivates 
a number of ambitious cluster surveys. Among these, surveys based on the Sunyaev-Zel'dovich 
(SZ) effect are particularly powerful for their ability to cleanly select clusters out to redshifts 
2 > 1. Now poised to begin surveying substantial areas of sky, dedicated interferometers, 
, bolometer cameras and the Planck satellite will soon produce large cluster catalogs that will 

04 ' provide a precise measure of the cosmic expansion rate over a range of redshifts and precipi- 

tate a new understanding of structure and galaxy formation. I review the science potential of 
, these surveys and examine some issues of SZ cluster catalog construction. 
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1 Introduction 
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Q ■ Detailed observations of the cosmic microwave background anisotropics , distant SNIa and 
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of the galaxy distribution ■ have driven the tremendous advance in recent years leading to 
^ I the development of a standard cosmological model. Ongoing research aims to test the model's 

coherence and to answer outstanding fundamental questions: What is the nature of dark matter 
and of the mysterious dark energy accelerating the present expansion of the universe? What is 
the physics of the early universe, in particular of the inflation epoch? How can we use detailed 



' cosmological observations to probe fundamental physics, such as the neutrino sector? And 

how do galaxies form and evolve in the cosmic web of large-scale structure. These exciting 
questions inspire large observational programs centered on the CMB, SNIa searches and wide- 
field multiband surveys. 

Galaxy clusters offer a unique avenue of attack on several of these key questions: 

1. Their abundance and evolution with redshift are highly sensitive to the statistics and 
growth rate of the density perturbations, which in turn depends on the cosmological pa- 
rameters; 

2. They are ideal tracers of the largest scale structures. With a mean separation of ~ 50 Mpc, 
they efficiently sample structures of wavelength ~ 100 Mpc and larger, such as the baryon 
acoustic peaks in the matter power spectrum (see Figure 

3. They provide a well-defined, quasi-closed environment for galaxy formation studies. In the 
cluster environment we directly observe the stellar, diffuse gas and dark matter components 
of the cosmic fluid; 



4. 



Combined X-ray and millimeter cluster observations (the Sunyaev-Zel'dovich effect; see 
below) permit distance measurements and hence the construction of a Hubble diagram with 
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Figure 1: Left: The effective volume as a function of wavenumber for the SDSS LRG sample (middle, red solid 
curve) and for two mass-limited cluster samples, one at M > 10^*Mq (upper, blue dashed line) and one at 
M > 2 X IO^^Mq (lower, blue solid line). In all cases the survey area is ^survey = 4000 deg^, corresponding to 
the actual survey area of the SDSS LRG sample that detected the baryon acoustic peak and to the proposed 
South Pole Telescope SZ survey, which targets these cluster mass limits. Right: The thermal SZ spectrum j^, as 
a function of dimensionless frequency x = hi^/kTcmb (Tcmb = 2.725 K). The zero crossing is at 217 GHz. 

which to measure the matter density and dark energy abundance and equation-of-state 
(hke SNIa distance measurements). 

2 Cluster Surveys 

Tapping the cosmological potential of galaxy clusters requires large, homogeneous cluster cat- 
alogs extending to redshifts greater than unity. The strongest constraints on dark energy, for 
example, will come from observed cluster evolution in the transition between matter and dark 
energy domination around z = 1 and beyond. Note that today there are only of order 10 clus- 
ters with spectroscopically confirmed redshifts z > l^^. The need for deeper and larger catalogs 
motivates a number of substantial observational efforts in different wavebands. 
We may search for clusters in a variety of ways: 

1. As overdensities in the galaxy population. This may either be done in redshift space, or in 
color space using multiband imaging which can provide estimated photometric redshifts 
over large areas of sky; 

2. By their gravitational lensing of background galaxy images. With a large imaging survey, 
we could search for clusters both via their lensing signal and in color space (as above); 

3. By the X-ray emission from their intracluster medium (ICM), which at keV temperatures 
emits copiously in X-rays. Cluster surveys and observations in the X-ray band with 
satellites like ROSAT, Chandra and XMM have significantly advanced our understanding 
of clusters and their evolution, and have played an important role in establishing the 
standard cosmological model ^. 

4. By the Sunyaev-Zel'dovich (SZ) effect, a distortion of the CMB frequency spectrum caused 
by a transfer of energy from the hot ICM electrons to CMB photons via Compton scat- 
tering. 



Each of these techniques has its own advantages and disadvantages. To appreciate them 
we should situate clusters in the standard theoretical framework that describes them as single, 
massive dark matter halos containing hot gas and galaxies with their subhalos. Dominated by 
dark matter and its gravitational evolution, clusters are fundamentally described by their halo 
mass and redshift. When building a cluster catalog, therefore, one should quantify its properties 
in terms of these basic parameters in order for proper comparison to theory. Specifically, we must 
accurately calculate the survey selection function and determine the observable-mass relation. 
Doing the latter is often referred to as calibrating the observable-mass relation, or for short, 
survey calibration. 

Large area redshift surveys are prohibitively expensive. Less expensive, imaging surveys 
rely on a consistent relation between the parent cluster halo and its member galaxy population; 
in particular, the relation between cluster mass and redshift, and galaxy colors and luminosity. 
Therein lies a potentially large uncertainty due to our limited understanding of galaxy formation 
and evolution (which these cluster surveys are designed to help us improve!). 

Lensing-based cluster surveys are appealing because they select clusters on their projected 
mass density. They will, however, suffer from the projection of mass along the line-of-sight, 
and their detection efficiency depends on the redshift distribution of the lensed galaxies. The 
latter is relatively easily controlled with photometric redshift information, while the im^ict of 
projection effects on survey selection functions and mass recovery requires further study . 

The properties of the hot ICM are closely linked to those of the cluster gravitational potential, 
implying that X-ray surveys select clusters based on their dark matter halos, similarly to lensing 
surveys. In addition. X-ray surveys do not suffer from projection effects, because cluster X-ray 
emission varies as the square of the gas density. For this same reason, however, the cluster 
X-ray emission is dominated by the core regions and hence physics that is difficult to model. 
Furthermore, X-ray surveys require large, expensive satellite missions. Not designed as survey 
instruments, Chandra and XMM-Newton will provide only relatively limited area surveys!^. 

Our focus here is on the fourth method based on the SZ effect, which possesses a number 
of advantageous properties for cluster surveying, as I detail below. As for X-ray observations, 
the SZ effect selects clusters based on their ICM properties and potential well; and like lensing 
surveys, SZ surveys tend to suffer from projection effects, since the SZ signal is proportional to 
the projected gas mass density. A number of planned SZ surveys have heightened anticipation as 
they prepare to observe large areas of sky over the next few years. These include ground-based 
instruments starting operations this year, as well as the Planck satellite scheduled for launch in 
2008 and which will cover the entire sky. 



3 The Sunyaev-Zel'dovich Effect 

The SZ signal arises from the energy transfer via Compton scattering from hot electrons in the 
ICM to CMB photons traversing the clusteJ^. Conserving photon number, the scattering diffuses 
photons from low to high energies and creates a unique frequency-dependent spectral distortion 
that is negative at low frequencies and positive at high frequencies. We quantify the distortion 
by the change in observed sky brightness towards the cluster relative to the unperturbed CMB 
seen from surrounding blank sky: 

= yj,y (1) 

where the Compton-y parameter specifies the amplitude of the distortion as an integral of the 
ICM pressure along the line-of-sight: 



(2) 



with T being the electron (ICM) temperature, m the electron mass, rie the electron number 
density and ctt the Thomson scattering cross section (we're in the non-relativistic limit of 
Compton scattering). The function gives the frequency dependence of the spectral distortion 
and is universal to all clusters at all redshifts (See Figure P). It is negative below 217 GHz, 
where the effect is a brightness decrement relative to the mean sky brightness (the CMB), and 
positive above, where the effect produces an excess brightness. 
We refer to the total SZ flux as 

Su = Yj, (3) 
where the integrated Compton parameter is 

4 SZ Cluster Surveys 

The SZ signal has a number of virtues when considered as a means of cluster surveying. Firstly, 
the surface brightness Ai^, is independent of distance and depends only on intrinsic cluster 
properties; this is also manifest by the fact that Y varies with the angular diameter distance, as 
opposed to the luminosity distance. Secondly, the spectral signature is unique, unlike any other 
astrophysical source, and it is universal to all clusters. 

Thirdly, the total SZ flux Si, is directly proportional to the total thermal energy of the ICM, 
a quantity fixed by energetics during cluster formation and independent of the final spatial 
and temperature structure of the gas. For this reason, we expect the SZ signal to tightly 
corre late to cluster virial mass, as indeed borne out in numerical simulations; for example, Motl 
et al.^^ simulated examined different simulations, including a variety of cooling and feedback 
mechanisms, and find only a ~ 5% scatter in y at a given mass. 

These characteristics of the SZ signal lead to one to expect that a SZ survey will select 
clusters based on their mass, and that the limiting mass will be nearly constant with redshift. 
This is a valuable property for studying any evolutionary trends, for we can compare similar 
mass objects at different epochs. 



5 SZ Catalog Construction 

In anticipation of the upcoming surveys, we have undertaken an extensive simulation effort to 
study survey selection functions and observational uncertainties, a nd I will report here two key 
results from this work to date. Details can be found in Melin et al^^ and Bartlett & Melin^J. 

We have developed a cluster detection algorithm based on a filter matche d to both the 
spectral (j^) and spatial (isothermal /3-model) signature of the cluster SZ signal The filter 
is designed to rapidly and optimally extract clusters of all sizes from the astrophysical fore- 
grounds. Rapidity is essential because we want to simulate a given survey many times in order 
to accurately quantify selection effects and observational uncertainties. 

Applying the detection routine to Monte Carlo simulations of the AMI, ACT, Planck and 
SPT surveys, we calculated the survey sensitivities, counts and completeness values shown in 
Figure [3 In this Figure, ay is the total noise through the matched filter generated by both 
instrumental noise and primary CMB anisotropy (for the standard WMAP cosmology), and 6c 
refers to the core radius of our simulated SZ profiles truncated at 10 core radii. The curves 
correspond to the different surveys as labeled. A cluster observed by SPT, for example, with 
angular core radius 6c and Y = ay {6c), i.e., on the solid, red line, would have a = 1. 

We see from this Figure that the surveys are not simply limited in flux; cluster selection 
depends on both fiux and angular size. This is particularly true for the high resolution ground- 
based surveys, while the lower resolution of Planck yields a nearly flux-limited catalog. In this 
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Table 1: Characteristics of some of the large SZ surveying instruments. We give the instrumental noise, counts 
detected at S/N > 5 with our matched filter and survey area for the 4 experiments that we have simulated in 
detail, namely AMI, ACT, SPT & Planck (see text). A list of web pages for these and other experiments is given 
in the reference section. Adapted from Melin et al. (2006) . Note added: The new WMAP3 results were 
published at the time of writing. The lower value of as favored by the new release can lower the predicted cluster 
counts by up to a factor ~ 2. We are evaluating the changes for each experiment in detail. 




Figure 2: The left-hand panei plots the matched filter noise cty as a function of filter scale 9c (core radius of a cluster 
matched to the filter) for different surveys, as labeled. The filter noise is generated by primary CMB anisotropy 
and instrumental noise. Clusters lying above the curve of a particular experiment have S/N > 1. Integrated 
source counts at S/N > 5 for each survey are shown in the right-hand panel, along with the simulation input 
counts (curve labeled "mass function"). We give the catalog completeness percentage (ratio of the experimental 
curve to the input mass function counts) in the inset. The important point is that the surveys are not flux limited, 
and are significantly incomplete even at 5 times their point source sensitivities (5 times the y-intercept of the 
left-hand curves). Adapted from Melin et al. (2006). Note added: The new WMAP3 results were published at 
the time of writing. The lower value of as favored by the new release can lower the predicted cluster counts by 
up to a factor ~ 2. We are evaluating the changes for each experiment in detail. 
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Figure 3: Photometric recovery in terms of integrated Compton Y parameter for SPT (left) and AMI (right). 
Compton Y values (in arcmin^ ) recovered by the matched filter are plotted against the input Y values taken from 
the simulation catalog. Each point represents a single cluster detected at S/N > 5. The red dashed curve gives the 
equality line. For SPT the characteristic scatter at fixed Ytruo is ~ 40%. Confusion with primary CMB anisotropy 
seriously compromises photometric recovery of the single frequency survey (chosen here as AMI). Adapted from 

Melin et al. (2006). 

Figur e w e only show the effects of instrumental noise and CMB anisotropies. In Bartlett &; 
Melin^^we examine the additional effects caused by extragalactic point sources. 

In the right-hand panel I show the integrated counts of clusters detected at S/N > 5. These 
are objects to be found in the left-hand panel above the corresponding curves, each displaced 
upward by a factor of 5. The counts roll away from the input counts (labeled mass function) at 
the faint end as the catalogs become incomplete. In the inset, we see the completeness defined 
as the ratio of detected to input clusters (ratio of one of the experimental curves to the mass 
function curve) . It is important to note that all of the surveys are significantly incomplete even 
at 5 times their point-source sensitivities, which are given by 5 times the y-intercept of the 
curves shown in the left-hand panel. This just illustrates again that the surveys are not simply 
flux limited. 

I now turn to another important, but often neglected issue - photometric recovery. In 
Figure ini I plot recovered Y versus input Y for clusters detected at > 5 in SPT (left- 

hand panel) and AMI (right-hand panel); the Planck case is similar to SPT. We see that when 
applied to SPT our algorithm gives an unbiased estimate of the true SZ flux, but with ~ 40% 
scatter (at fixed i^rue)) which is much larger than the intrinsic scatter predicted by numerical 
simulations^!. The photometry of the single frequency AMI survey, on the other hand, is 
very seriously compromised; the distribution at fixed Itrue is in fact bimodal. We have traced 
this effect to confusion with primary CMB anisotropy that is unavoidable with single frequency 
observations. The conclusion is that such single frequency surveys must follow-up their clusters 
at a second frequency in order to eliminate the confusion. 



6 Conclusions 

Future galaxy cluster surveys will provide a wealth of information on dark energy, dark matter 
and structure and galaxy formation. Among these surveys, those based on the SZ signal will 
profit from its intrinsic ability to find clusters at high redshift and its expected tight correlation 
to halo mass. Over the next 5 years, these surveys will provide large, well-defined (in terms 
of mass) catalogs containing hundreds to thousands of clusters at redshifts beyond unity or, in 



other words, multiplying by 10-100 the number of known clusters at these redshift^. This will 
give us a new view of galaxy formation in dense environments and a measure of the expansion 
rate at this crucial epoch marking the transition between matter and dark energy domination. 

Many of the surveys will begin this year (2006), and a number of important surveying 
issues require further study. These include accurate evaluation of survey selection functions 
and observational errors. I have shown two results from our studies based on a matched filter 
detection algorithm and simulations of different surveys. We find that high resolution ground- 
based surveys select clusters not simply on flux, but on a combination of flux and angular size, 
and that this must be properly accounted for when modeling cluster counts. We also find that 
observational scatter on measured cluster flux Y (e.g., ~ 40% in the case of SPT) largely exceeds 
the intrinsic scatter predicted by numerical simulation^^. Furthermore, confusion with primary 
CMB anisotropy severely compromises photometry in single frequency observations; follow-up 
at a second frequency (or in X-rays) will therefore be necessary for these surveys. 
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• ACBAR http : //cosmology .berkeley. edu/group/swlh/acbaj/l 

• ACT jhttp : //www . hep . upenn . edu/^eingelica/ act /act . htmlj 

• AMI [http : / / www . mrao . cam . ac . uk/telescopes/ami/index . htinl| 

• AMiBA [http : / / www . asiaa . sinica . edu . tw/ amibaj 
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• SZ A jhttp : //astro . uchicago . edu/ sze| 

• Olimpo fhttp: //oberon.romal . infn. it/| 
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